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Abstract 

Exchange of unparticle stuff of dimension du with FCNC interaction can induce M 12,u and r 12,u 
causing meson and anti-meson mixing with the relation T 12,u /M 12 ' u = 2tan(7rc£^). We show that 
this type of unparticle contribution can provide the much needed large T 12 to explain the recently 
observed anomalously large dimuon asymmetry in B s — B s system reported by DO collaboration. 
The same interaction can also accommodate large mixing induced CP violation in B s — > J/ip(fi 
indicated by CDF and DO data. Experimental data can provide constraints on the unparticle 
dimension and scale. 
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Recently the DO Collaboration has reported evidence for an anomalously large CP viola- 
tion in the like-sign dimuon charge asymmetry [1] which is attributed to semileptonic decays 
of b hadrons defined by 

A b = ^ ++ ~ (1) 

sl " N++ + N-- ' K } 

where N^~ + (N b ) is the number of events with two b hadrons (b b) decaying semileptonically 
into fi+fj+X {n~n-X). The DO result A b sl = -(9.57 ± 2.51 ± 1.46) x 10~ 3 with the first 
and the second errors being statistical and systematic ones, is 3.2a away from the standard 
model (SM) prediction of -0.2 x 10" 3 [2]. 

A^ is related to the "wrong-charge" asymmetries aff in and B s decays 

q = T(B q -» fi+X) - T(B q -» f£X) 
Gsl T(B q -)■ fi+X) + T(B q -> fi-X) ' 1 } 

With known values for mixing parameters of Bd and B s systems, one obtains 

A b sl = (0.506 ± 0.043)aJ + (0.494 ± 0.043)a s S ! . (3) 



The above leads to LU af, = —0.0146 ± 0.0075, after insert the known value [lj, \A\ 

n ft 

0.0047 ± 0.0046. Combining direct information on a s sl from CDF [5j and DO [6j, one can 



mbin 

extract an average value of 



4 < 



K)ave~ -(12.7 ± 5.0) xlO' 3 . (4) 



The "wrong- charge" asymmetry for B s — B s mixing is determined by the matrix elements 
M s 12 and T] 2 in the Hamiltonian, 

Im(Mj 2 n 2 *) _\T] 2 \ . 

G sl" | M 12|2 + | r i2|2 /4 W |Mi 2 | Sm0S - 

In the above, for the last step of approximation, we have neglected small corrections of order 

|r?| 2 /|M s 12 | 2 . 

The SM prediction of |r] 2 ' SM |/|M s 12 > 5M | is (4.97 ±0.94) x 10~ 3 and the phase S is small 
leading to a value (2.1 ± 0.6) x 10~ 5 for a|j [2[. This is 2.5cr away from (a^ave, and the 
predicted central value is more than two orders of magnitude smaller. If the DO result is 
confirmed, it is an indication of new physics beyond the SM. Implications of this result have 
been studied recently by several groups 
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Different values of M] and T] than those predicted by SM are needed. In many models 
beyond the SM, there are new FCNC interactions which can easily induce large modification 
to M} 2 , but not to T] 2 , such as Z 1 and SUSY models. If modifications to T] 2 is negligible, 



i.e. r 
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(6) 



(7) 



In the above, M^ 2 ' NP indicates the contribution from beyond the SM new physics. We have 
neglected small SM contribution to <p s and adopted the phase convention that T 12 is real. 

Since the SM contribution to the phase <p s is small, it is convenient to write the ex- 
pression explicitly in terms of the new physics contribution defined by M} 2,NP /\Ml 2,SM \ = 
R exp(—i(f) NP ), one then has [17 ] 



U2,SAT| 



i?sin 



lNP 



| Ml 2 ' SM I 1 + 2R cos <P NP + R 2 
\A S \ = v 7 ! + 2R cos <f) NP + R 2 . 



(8) 



The SM prediction for AM S w 2\M] 2 \ agrees with data well, |A S | is only allowed to vary 
in a limited region, 0.92 ±0.32 fixed by the experimental value [3J] AM S = (17.77 ±0.12)ps~ 1 
and the SM prediction [2] (19.30 ± 6.74)ps _1 . In order to reproduce the DO result, naively it 
seems that sin <p s would have to exceed the physical range | sin <p s \ < 1 to reach the central 
value from data. A large beyond SM new physics contribution Y l2,NP to V 12 is therefore 

n n " 

called for. Several attempts have been made in this direction |9|, [16|, 



17 



30]. 



In this work we show that unparticle FCNC interaction can provide the much needed 
large r] 2 ' n to explain the observed anomalously large dimuon asymmetry in B s — B s system 
by DO collaboration. Using experimental data, one can constrain the dimension du and the 
unparticle scale. 

Several years ago Georgi proposed an interesting idea to describe possible scale invariant 



effect at low energies by unparticles 



341 ] . An unparticle operator Ou of dimension du result- 



ing from some scale invariant theory at high energy may interact with SM particles at low 
energy in the following form 



U U SM U U- 



(9) 



where Osm is composed of SM fields. 

If the Osm is a bi-quark operator with different quarks, exchange of unparticle stuff can 
induce meson and antimeson mixing [35]. Unlike usual tree level contributions to meson 
oscillations from heavy particle exchange which produces a small T 12 , the unparticle may 
have sizeable contributions to both M 12,u and T l2 ' u due to unparticle fractional dimension du 



leading to a phase factor 34j, |36) (— \) du ~ 2 = e~ mdu in the propagator which is proportional 
to M 12 — iT 12 /2. This leads to an interesting relation y 12 > u /M 12 ' u = 2tan(irdu)- Depending 
on the dimension du, T 12,u can even be larger than M 12,u . 

Contributions to M l2,u and T l2,u for a class of operators Osm with dimension less than or 



equa 
Ref. 



to 4 interacting with unparticles, respecting SM symmetries 37|, have been studied in 



Here we take a scalar unparticle Ou interaction with quarks, respecting SM gauge 



symmetry, of the following forms for illustration, 

zAgA^4 T X<W w , fX^A^QhM^Ou , (10) 

where d l R and Qf[ = {u\, d l L ) are the right-handed down-quarks and left-handed quarks, 
respectively. 

After using equation of motion for quarks, one can write the above tree level FCNC 
interaction as 

C = ^ d ^m^ R ~ L ) + *g(^ - R )]djO u • (11) 

Let j be the heavier and i be the lighter quarks, the term proportional to m;Jmj can be 
neglected. 

Note that the above Lagrangian, in general, can induce B s — > Ou transition making a 

large contribution to B s invisible decay rate for i = s and j = b. If there is no similar 

interaction for i = d, this will also induce a significant lifetime difference between and 

B s . These problems can be solved if the couplings \® b and A^, are equal, A£ = A^, such that 

the above Lagrangian is proportional to sbOu which cannot induce B s — > Ou transition. We 

will work with this condition. 

Exchange of unparticle in the t and s channels as shown in Fig. 1, we obtain the effective 
tU 



Hamiltonian H^ ff responsible for meson-antimeson mixing 
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H eff = o A <\ ^ u'^^-A (l \xTir + a L I ) ( d ^S L + ~ X % R ) d i)h) 
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FIG. 1: The t and s channel contributions to meson-antimeson oscillation. 



Here A d , 



(167r 5 / 2 /(2vr) 2 ^)r(^ + l/2)/{Y{d u - l)T(2d u )). 



We have used 



{iA du /2 sm( y 7idu))e~' ldu7T (l/(p 2 — n 2 ) 2 ~ du ) for the scalar unparticle propagator [39|. One 
can easily identify the factor e~ tndu in the propagator producing a non-zero r 12 '" mentioned 
earlier. 



In the above we have followed Ref. [391] to use a modified unparticle propagator to take 
into account the breaking of conformal invariance at a scale \x lower than An by unparticle 
interaction with SM particles. Whenp 2 becomes smaller than /i 2 , the unparticle contribution 
to T 12 vanishes. The scale \x is not known, since we rely on the unparticle effect to generate 
T^ 2 , we will allow it to be lower than B s mass. One should not allow /x to be too small to 
allow b — > sOu process to happen which will contribute too large an invisible decay rate for 
b sOu which changes B s lifetime dramatically. For this reason we assume that \x is in the 
range of > // > ra^ — ttik such that b —> sOu is forbidden. 

In B s — B s system, mesons are made of a i = s light quark and a j = b heavy. In the 



heavy quark limit, one has s = t 



mi 



m 2 B . With this approximation and theoretical 



matrix elements for the relevant operators, we have 



M s 12 '" - i-T 12 ' u 
2 



.4 



du 
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x ( (Af b ) 2 + (A«) 2 - 6((A£) 2 - (A«) 



fi 
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12m, 
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ZA sb A sb 



(13) 



where Bb s is the bag factor which is equal to 1 in the vacuum saturation and factorization 
approximation. With the condition A£ = \® b , the factor in the second line of the above 
equation becomes 4(A^) 2 . 

It is clear from Eq. (Tl3l) that because of the phase factor e~ mdu in the propagator, there 
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is a relation between M] 2 ' u and T 12,u that 

^- = 2tan(nd u ). (14) 

There may be a sizeable contribution to T 12,u at tree level which is not possible for usual 
tree level heavy particle exchange. For dy equal to half integers, there is no contribution 
to M s 12,u , but there is for T 12,s . It has been shown 38j] that a large class of unparticle tree 



level interaction contributions to meson and antimeson mixing can be written in the form 
of Eq. (fT51) . 

One can group unknown parameters in such a way that 

M} 2 ' u = |M, 12 ' M |e-^ , r s 12 -" = 2|M s 12 '"|e-^ tan(vrrf w ) , (15) 

where 

|Mf>-^ = \ Adu U Ji±-B Bs cot(d M 7r)|e-^ . (16) 

Here we have identified the phase of A^ to be —0^/2 and defined the reduced unparticle 
scale A u = A u (l — /i 2 /m 2 Bs ) l l du ~ 1 l 2 \\® b \~ l / du . ~k u is an effective indicator of the scale of 
unparticle physics. This quantity actually contains three unknown quantities, Ay, |A^J and 
fi. Since we do not have good ideas about their individual sizes, we will use the reduced 
scale for discussions later. In general A^, is complex, therefore <pu can take an arbitrary 
value between to 2n. 

Combining the SM contribution and neglect small phase there, we have 

s = \T 12 ' SM \ flsinQfo) _ 2taxi(dun)R S m(<ki) 

fl si | M i2,SM| ! + 2 J Rcos(</» w ) + R 2 1 + 2Rooa{<ki) + R 2 ' 

|A S | = v /l + 2 J Rcos(0 w ) + R 2 . (17) 
The first term in the expression for is similar to the contribution without modification to 



T^ 2 analyzed in Ref . 17| . The second term is new which comes from unparticle modification 
to I?. 

There are two constraints with three parameters, R, (pu and y = tan(c/^7r). It is not 
possible to determine all the theoretical parameters. However, one can express two of the 
theoretical parameters as a function of the experimental measurable quantities and another 
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theoretical parameter. We find 

COS(0 W ) = - 



2R 

The "±" sign in the expression for y comes from the need of taking a square root for sin(0^) 
from knowing cos(0w). 

Since the physical range of cos(0^) is limited to be from -1 to 1, from the first equation 
about, one can obtain a constrain for the allowed ranges for R. We show this in Fig. 2 
for several values of |A S |. In the plot we have used the central values for the relevant SM 
quantities. 

With the allowed range for R known, for given and |A S |, y can be determined as a 
function of R. For example, with central values for and |A S |, for R = 0.5, 1.0, 1.5, we have 
y = 0.008,0.009,0.012 and y = -0.003,-0.004,-0.007 with "-" and"+" signs in eq. 18, 
respectively, and |A S | can be accommodated with reasonable values for other parameters. 
y can take both positive or negative signs. To have further information on d u , one should 
also consider dy = n + (arctan y)/n with n an arbitrary integer. However, if one considers 
the relation between the reduced scale k u and dy, one would prefer not too large a du- If 
one limits du between 1 and 2, when y takes a positive value, one should take n — 1. When 
y is negative, it is reasonable to take n = 2. 

The unparticle interaction we have introduced can also have effects on other flavor chang- 
ing processes which may exclude some of the parameter space favored by producing the 
value obtained by DO. For AS* = 1 processes, such as B S ( i decays with SM final states and 
unparticle as intermediate state, there is the need to introduce new unparticle interaction 
with other SM particles, one can adjust the new couplings to be sufficiently small to satisfy 
bounds from experimental data. This type of interaction will not affect the analysis carried 
out earlier and we will not discuss them further. There are, however, experimental observ- 
ables which depend on the quantities we have discussed. For example mixing induced CP 
violation in B s decays can be affected by modifications to M s 12 , in particular, experimental 
data on mixing induced CP violation in B s — > J/^ip can further constraint the interaction 
discussed. 

Recently analysis from CDF and DO data show that there are indications of mixing 
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0.0 0.5 1.0 1.5 2.0 

R 



FIG. 2: R vs. cos(^). The solid and dashed curves are for |A S | and sinc/> s taking their central 
and one a upper and lower values, respectively. The shaded areas are the one a overlap regions. 
(|A S | = 0.6 (a), | A, | = 0.92 (b), |A S | = 1.24 (c); sin^ = -0.470 (A), sinc^ = -0.704 (B)), 
sin^ s = -0.908 (C)). 




FIG. 3: An vs. du for fixed R = 0.46; 0.73; 1.78; 2.00, in order from up to down (magenta (a); 
green (b); blue (c); red (d)), respectively. 



induced CP violation in B s — >■ J/ip<p decay with [40] 



ft. 



(o.39±8:i5) u (1.181S1S) 



(19) 



where the second set in the last line is just the complement of ir/2. This reflects the ambiguity 



in the determination of fi 3 J^ . 



If there is no new interaction in the decay amplitude, —2(3$^^ is the phase of M s 12 . 



The SM prediction for is given by |2|, /3 s w(SJVf) = Arg(-^gl) « 0.019 ±0.001. The 

central value of ft J J^^ SM > i s much smaller than that from CDF and DO and can be neglected. 
If confirmed, this also requires new physics beyond SM. Even though we do not introduce 
new unparticle interactions in the decay amplitude, M s 12 is modified which can affect mixing 
induced CP violation in B s — > J/tpcj). We have —2fi J J^"^ ~ 4> s . Here <p s is defined in Eq. (JTj). 
The two values for /3/ give approximately the same sin(0 s ) with the central values given 
by -0.704 and the one a range of -0.908 ~ -0.470. 
Translating into R and <pu defined earlier, we have 

. -Rsin(<f) U ) 

smd) s = — . 20 

v /l + 2 J Rcos(0 w ) + J R 2 

It is clear that a non-zero <pu can modify the SM prediction for /?/ . 

Using the above one can further constraint the allowed regions for R and cos(0^). From 
the relation between (p s and 4>Ui one finds sin(0^) must be positive. The allowed regions for 
cos(0^) are shown in Fig. 2. The shaded regions are the one a allowed common regions. 

With three unknowns and three input experimental data points, we can completely ex- 
press R, cos(0^) and y as functions of |A S |, (fi s and a s sl . We have 



R= Vl±2|A s |cos0 s + |A S | 2 , 

(A \ ~ lT |A s |cos0 s 

cos((pu) = , , (21) 

v/l±2|A s |cos0 s + |A S | 2 

v = -l L + sl1 s 



2 V|M S 12 ' 5M | sin0 

We now use the above equations to determine central values and their one a errors for 
R, cos(0^) and y. Since the error correlations between the three experimental observable 
quantities are not know, a complete analysis is not possible. To have some idea about the 
ranges for the theoretical parameters, we will assume that they are independent quantities 
satisfying Gaussian distribution and take the average one a lower and upper bounds as their 
one a errors for illustration. The results are listed in Table 1. For a given set of input data, 
there are two solutions for R and cos(0^), but the solution for y is single valued. We see 
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R 


cos{4) U ) 


y 


1.776 ±0.185 


-0.931 ± 0.055 


0.011 ±0.004 


0.735 ± 0.276 


-0.472 ± 0.208 



TABLE I: The central values and one a bounds of R, cos(4>u) and y. 



that y is restricted to be positive now. The large dimuon asymmetry in B s — B s mixing and 
large mixing induced CP violation in B s — > J/ip(j) can be simultaneously explained. 

Finally, we briefly comment on the implications of constraints on the effective unparticle 
scale An. Using the results obtained on R, we can work backward to obtain information 
about the reduced unparticle scale, A^. For a given R, An can be expressed as a function 
of du- We show the results in Fig. 3. In plotting Fig. 3, we have used B Bs = 1 and 



0.260 GeV 



4l| . We can see, as expected for larger du-, a lower scale is needed which 
may have already reached by current colliders. Since no evidence of unpaticle effects have 
shown up explicitly, one should set a lower limit for that. For a lower bound of A^ larger 
than 1 TeV, we see that du is restricted to be lower than 1.5. 
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In summary we have shown that unparticle induced contribution to B s — B s mixing 
can provide the much needed large V 12 to explain the recently observed anomalously large 
dimuon asymmetry by DO collaboration. The same interaction can also accommodate large 
mixing induced CP violation in B s — > J/ip<p from CDF and DO data. Experimental data 
also provide constraints on unparticle dimension and scale parameters. 
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Note Added 

n 

While we are in the final stage of finishing our work, a paper |42| by J. P. Lee appeared 
where unparticle effects on B s — B s was also discussed in the context of large dimuon 
asymmetry reported by DO. We do not agree Lee's interpretation for the part proportional 
to sm(7idu) generated due to the e~ tndu factor in the propagator. It contributes to T] 2 , but 
not M] 2 . Our emphasis is on how the unparticle effects modification for T] 2 plays a crucial 
role in explaining the data reported by DO. In general a CP violating phase, in our notation 
(fin, can also exist. This phase is not considered in the analysis by Lee. 
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